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We report a new measurement of the decay f2 



£1 7T T 7T 



with 24 events, yielding a combined branching ratio (3.74_ggg) x 10 



with 76 events and a first observation of the decay 
+0 - 67 ^ - ln ~ 4 . This represents a factor 25 



increase in statistics over the best previous measurement. No evidence is seen for CP violation, with B(tt 



— > H 7T + 7T 



4.04 



+0.83 
-0.71 



x 10" 4 and 



3.15 

— *o 



+ 1.12 
-0.89 



x 10 . Contrary to theoretical expectation, we see little evidence 



i for the decays Q ^i® 30 Tr and Q + — > ^ 1530 7r + and place a 90% C.L. upper limit on the combined branching ratio 
B(fi-(n+) -> H^o^&o)^) < 7-0 x lO" 5 . " 
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Although the Q~ hyperon was discovered in 1964, many 
of its decay modes are still poorly known experimen- 
tally. One of these is Q~ —> E~ir + ir~. The f2~ is the 
only ground-state hyperon massive enough to have such a 
AS* = 1, three-body, nonleptonic weak decay. This decay 
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has long been assumed by theorists to proceed predomi- 
nantly via the £*5 30 intermediate state [TJ[2]. Based on 



four observed ft — > S 7t + tt events, Bourquin et al. es- 
tablished the best previous measurement of the branching 
ratio: B{Qr -> E-jt+tt-) = 4.3±?;| x 10~ 4 [3J. However, 
they were unable to distinguish the nonresonant chan- 
nel from the resonant 0~ — > H*53 7r~ — » S~7r + 7r - de- 
cay. In interpreting this observation, the Particle Data 
Group (PDG) adopted the assumption by Finjord and 
Gaillard [2] that the S~7r + 7r~ final state in Q~ decay arises 
entirely from the S^g 30 , hence that B(fl~ —> S~7r + 7r~) = 

|B(n--> 

mate B(Q' 



7r ) by isospin symmetry, giving the esti- 



"*0 
"1530 
- v ^7*0 - 
^ "1530~ 



) = 6.4 



+ 5.1 



10~ 4 [I] . Yet, until 



now, that assumption has not been tested by experiment. 

The HyperCP experiment (Fermilab E871) took data 
in 1997 and 1999. The apparatus, depicted schematically 
in Fig. [I] is described in detail elsewhere [5]. In brief, 
an 800GeV/c proton beam struck a target located at the 
entrance of a collimating channel within the "hyperon" 
dipole magnet. The selected momentum ranged from 110 
to 240GeV/c, with an average of rj 160 GeV/c and r.m.s. 
spread of s» 25 GeV/c. A 13-meter-long evacuated decay 
volume started at the exit of this hyperon channel and 
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Figure 1: Plan view of the HyperCP apparatus. The typical 
topology of the signal decay fi — — > H~7r+7r~ — > (Att~ )7r+7r _ — > 
(ptt~tt~)it + it~ is shown. 

was viewed by a charged-particle tracking spectrometer 
comprising four multiwire proportional chambers preced- 
ing and five following a pair of dipole analyzing magnets. 
The polarities of the magnets were such that particles hav- 
ing the same charge as the beam or opposite to it were 
deflected to the left or right, respectively. The tracking 
chambers were followed by trigger hodoscopes in each arm 
of the spectrometer. Also used in the trigger was a hadron 
calorimeter in the "opposite-sign" arm. For the data sam- 
ple discussed here, the trigger requirement was at least 
60 GeV deposited energy in the calorimeter and hits in 
the hodoscopes consistent with the passage of at least two 
charged particles of opposite polarities. Antihyperon pro- 
cesses were studied using data recorded under identical 
trigger conditions (referred to below as "positive" data) to 
those for hyperons by reversing the magnetic-field direc- 
tions of both the hyperon and analyzing magnets; singles 
rates in the spectrometer were approximately equalized by 
using a target one-third as long for positive data. 

The experiment recorded ^30,000 5 GB data tapes con- 
taining some 2.3 x 10 11 events. These raw tapes were pro- 
cessed [6] to produce output data sets satisfying various 
topological or particle-ID requirements depending upon 
their intended use. The selection cuts for the analysis 
presented here were initially devised from Monte Carlo 
(MC)-simulated data samples for signal (f2 _ — > E~ir + ir~, 
S~ — > Att~ , A — » p7r~) and normalizing — > AK~ , 
A — > pir~ , K~ — > 7r~7r + 7r~) modes, prior to examination 
of the data samples. The cuts used in this analysis re- 
quired: (a) at least three negative and two positive tracks 
(or, for the positive data sample, three positive and two 
negative tracks) [7], (b) fi, S, A, and K decay vertices lo- 
cated inside the decay volume, (c) a topology consistent 
with the intended decay, (d) a combined 5-track momen- 
tum between 135 and 220 GeV/c, (e) reconstructed Q, S, 
A, and K masses within ±3 standard deviations (a) of 



their nominal values [3], and (f) that the reconstructed SI 
track extrapolate to the target and pass through the aper- 
ture of the collimator. By Monte Carlo simulation we find 
the overall efficiencies of these cuts to be 50% for the sig- 
nal mode and 43% for the normalizing mode [8] . The MC 
simulation was carefully checked against data and success- 
fully reproduced mass resolutions and particle spatial and 
momentum distributions. 

In each event, the opposite-sign track with the highest 
momentum was assumed to be the proton. MC simula- 
tion showed this assumption to be correct for more than 
99% of events. This assumption also gave a reconstructed 
invariant mass closest to the PDG value, when trying 
all possible track combinations, for 100% of data events 
with five charged tracks and for 98% of events with more 
than five charged tracks (which constituted ~ 35% of the 
data sample). In determining which same-sign track to 
associate with the proton to form a A, that pion resulting 
in the best A mass (i.e., closest to the PDG 2008 average 
value of 1,115.68 MeV/c 2 j3]) was chosen, and similarly 
for the H. For normalizing-mode events with more than 5 
tracks, after reconstruction of the A, the three remaining 
tracks of the appropriate polarities giving the best K mass 
were used. 

We normalized the candidate events using the — » 
AK~ , A — > pn~ , K~ — > 7r~7r + 7r~ decay since its topol- 
ogy is similar to that of the signal mode. The MC was 
used to calculate the HyperCP apparatus acceptance cor- 
rections for the normalizing and signal modes. The signal 
mode was simulated with a uniformly populated phase- 
space generator. 

Results based on the larger, 1999 HyperCP data sam- 
ple are reported here. The observed negative (positive) 
five-track mass distributions of the signal and normaliz- 
ing decays are shown in Fig. [2] (Fig. [3]); we observe 78 
(24) events within ±3<r of the MC mass resolution (a = 
2.0 MeV/c 2 ). The unbinned log-likelihood fits shown in the 
figures employ a Gaussian for the mass peak, with r.m.s. 
width (2.0MeV/c 2 ) determined by fitting the MC signal 
and normalization samples, plus a constant representing 
the background. From the fits, we estimate 2.3 ± 0.5 
(0.1 ±0.1) background events under the peak, within ±3<r, 
for the negative (positive) signal mode. For the normaliz- 
ing mode, 375 (156) events are observed within ±3cr of the 
MC mass resolution including an estimated background of 
0.5 ±0.3 (0.4 ±0.2) events. 

The signal events presumably represent a combination 
of direct Q~ — > 'Er tt + h~ decays with branching ratio Bdir 
and those proceeding via the resonant decay mode iY~ — > 
^i530 7r ~ with branching ratio B res , with the subsequent 
decay — > E~ir + . Because of the short lifetime of the 
2^530, the Sj53 and S~ vertices are indistinguishable. The 
main difference between the resonant and direct modes 
that we have identified in MC studies is the respective 
distributions in S _ 7r + invariant mass: as shown in Fig.[4j 
the resonance-mode mass distribution is sharply peaked at 
the resonance mass, while the direct one is more broadly 
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Figure 2: Observed-event histograms for (left) SI — > 3 7r+7r signal mode and (right) normalizing mode f! — > AK where the A decays 
into pn~ and the K~ decays into three charged pions; dashed curves show Gaussian-plus-constant fits to the data as described in text. 
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Figure 3: Observed-event histograms for (left) f2+ — > S+7T+7T signal mode and (right) normalizing mode H+ — > KK + where the A decays 
into p-7T+ and the K + decays into three charged pions; dashed curves show Gaussian-plus-constant fits to the data as described in text. 
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peaked at rj 1.52GeV/c 2 . We decompose the N sig signal 
events into numbers of direct (Ndir) and resonance {N res ) 
events according to 



N, 



dir 



Nn-->E-7T+7T- = fdir X N, 



N res = N n - 



sig , 



f, 



r<ia ^ sip ■ 



(i) 

(2) 



Table 1: Direct 3-body and resonance- mode fractions of observed 
Q~ — > H~7r+7r~ and f2+ — > signals, based on generalized 

likelihood fits to data. 



Decay 


fi- 


n+ 


fdir 


0.851±0.128 


0.973±0.245 


fres 


0.149±0.086 


0.027±0.144 



To determine the resonance {fres) and direct {fdir) frac- 
tions we fit the signal-event S~7r + mass distributions with 
a linear combination of the functional forms shown in 
Fig. [4] Results are summarized in Table [T] the fit functions 
together with corresponding experimental S~7r + invariant 
mass distributions are plotted in Fig. [5j Due to the small 
numbers of observed events, the fits were performed us- 
ing the unbinned generalized log- likelihood technique [9]. 
The goodness of fit was evaluated using the x 2 and the 
Kolmogorov test and was acceptable in all cases [8 ]. No 
significant resonance component was observed. 
We calculate 



B{Q- 



"1530 7r ) — fr 



N ■ A 

ly sig -^norm 



Bn- 



>AA'- 



N A 

ly norm 
■ B K -^k+ 
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7T T 7r — 7T — 



B^o 

^1530 



— >At> 



(3) 



and also obtain signal-mode branching fractions for the 
direct process 



Bdir{tt~ 



5 7T+7T ) = fdirjj 

Bq-^ak- ' B K 



N sig A nor m 



X ■ 



►Att" 



(4) 



where N is the number of events and A is the accep- 
tance for the indicated mode as derived from MC simu- 
lation ^1 The known branching ratios entering into Eqs. [3] 
and [ZJare B(Qr AK~) = (67.8 ± 0.7)%, B(E~ 
Att-) = (99.887 ± 0.035)%, £(S*° 30 -> E-JT+) = 2/3, and 
B(K~ -> ir-TT+n-) = (5.59 ±0.04)% T. The numbers 
of signal- and normalizing-mode events are obtained from 
the Gaussian-plus-constant fits of Figs. [2] and |3j by sub- 
tracting the fitted number of background events from the 
number of events observed within ±3<r of the mass peak. 
Table [2] summarizes the relevant quantities and results. 
Since we cannot identify which signal-candidate events un- 
der the mass peak are background and which are signal, 
we performed a study in which randomly chosen candidate 
events were excluded to see how much the results changed; 
as expected given the large signal-to-background ratio, the 
changes were negligibly small. 

Systematic uncertainties arise from uncertainties in the 
known branching ratios of secondary or normalizing decays 
as reported in [3] , variations and uncertainties in the appa- 
ratus magnetic fields and in the target and beam locations 



13 Note that the presence of B(Hjg 30 - 
of Eq. |3] takes into account the mode 
analysis is insensitive. 
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) in the denominator 
• 5°7r , to which this 



Table 2: Observed numbers of events, estimated background, accep- 
tances, and measured branching ratios; errors are statistical only. 



Quantity 



fi- 



N° bs 

norm 
ivrbkg 
norm 

N° bs 

sig 
lybkg 

sig 

A 

-fir/.orrn. 



Adv 
A 

si-re 

E 
B, 



dir 



375 
0.5±0.3 

78 
2.3±0.5 
2.81 x 10~ 4 
4.88 x 10~ 3 
1.20 x 10" 2 
(0.40 ± 0.24) x 10" 
(3.75 ± 0.74) x 1Q- 



156 
0.4±0.2 

24 
0.1 ±0.1 
2.78 x 10~ 4 
4.92 x 10~ 3 
1.20 x 10~ 2 
(0.06 ± 0.29) xl0~ 4 
(3.20 ± 1.07) xlQ- 4 



during the run, and uncertainties in relative geometric ac- 
ceptance between signal and normalizing modes due to 
uncertainties in the £1 momentum spectrum and the f2, 2, 
and A decay parameters [J]. Although the results of our 
analysis are dominated by statistical uncertainty, system- 
atic errors were conservatively estimated using the MC, by 
individually varying the respective parameters by ±(>1<t), 
and recalculating the branching ratios with the modified 
parameter value (details may be found in Ref. [5])- These 
error estimates are summarized in Table [3j where (to be 
conservative) we have symmetrized some asymmetric er- 
rors using the worse case. Because this analysis is nor- 
malized using a decay mode with five charged particles 
obtained using the same trigger and offline data set as the 
signal modes, systematic errors due to such effects as trig- 
ger efficiencies tend to cancel and are therefore neglected. 

Contrary to theoretical expectation [3J [TU] , we see little 
evidence for the resonant decay. To derive an upper limit 



Table 3: Systematic errors (in %) by source and data sample. 





(6B/B) dir 


{8B/B)res 


Source of Error 


n- 


0+ 


fi- 


n+ 


Branching-ratio and 


4.20 


4.20 


3.54 


3.54 


lifetime errors 










MC momentum spectrum 


3.78 


3.78 


6.79 


6.79 


Beam targeting 


3.38 


3.38 


3.54 


3.54 


Magnetic fields 


1.69 


1.69 


1.84 


1.84 


Decay parameters 


1.94 


1.78 


2.08 


1.93 


MC statistics 


1.30 


1.30 


1.30 


1.30 


Total (%) 


7.19 


7.15 


8.97 


8.94 
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Figure 4: Mass distribution of pn + n it for (left) resonance-mode and (right) direct decays from MC simulation; curves are analytical 
functions determined by smoothing these histograms. 
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Figure 5: Mass distributions of (left) pir + n n and (right) pn 7r + 7r + , together with fits (solid curves) to sums of resonance-mode (dashed 
histograms) and direct (dotted histograms) functional forms; dashed and dotted histograms are each normalized to the total number of events. 
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on the resonant contribution including both statistical and 
systematic uncertainties, we performed a Monte Carlo sim- 
ulation of a large sample of hypothetical experiments that 
took into account the probability distribution from the 
unbinned maximum-likelihood fit as well as the uncer- 
tainty in the branching-ratio normalizing factor (treated 
as Gaussian-distributed). We find 

B(n- -> s;g 30 7r-) < 7.4 x 10- 5 , 

B(Ti+ -> S^ 30 7r+) < 6.2 x 1(T 5 , and 
B(fi* -> Zf 530 (Zf 530 )nT) < 7.0 x 10- 5 

at 90% confidence level (C.L.). 

The total branching ratio for decaying to S~7r + 7r~ 
was calculated as a sum of Bdir and B res x E(3ig 3 Q — > 
S~7r + ). Monte Carlo simulation was employed (as de- 
scribed above) to combine statistical and systematic er- 
rors taking correlations into account. Median values to- 
gether with the ±34.1-percentile deviations for the re- 
sulting asymmetric distributions of branching ratios are 
B(fi~ -> S-tt+tt^) = 4.04±{j$ x 10~ 4 and B(U+ -> 
E+iT+n-) = 3.15±J'jg x 10~ 4 . We see no evidence for CP 
violation: particle and antiparticle branching ratios are 
consistent with each other. Defining the CP asymmetry 

= B(n- -» E-TT+TT-) - B(U+ -» S + 7T+7T-) 

= B(n- -> S-7T+7T-) + £(H+ -> S+7T+7T-) ' ^ ' 

we find ^4 = 0.12 ± 0.20. Combining results for par- 
ticles and antiparticles, we obtain an average branch- 
ing ratio B(nT -> E^tt+tt-) = (3.74±gjJ) x 10" 4 . 
In contrast to the expected |, this implies B(D, T —> 
Et° 30 {El° 30 )TT^)/B(n^ -> S^tt+tt-) < 0.20 at 90% C.L. 

In conclusion, using a 25-times-larger sample than that 
of the previous measurement [3], an improvement in pre- 
cision has been achieved in determining the branching ra- 
tio for SI - — > S~7r + 7r~. Our observations of the f2 _ and 
f2 + modes agree with CP conservation and the final com- 
bined branching ratio is (3. 74^'^) x 10~ 4 . Contrary to 
expectations, we see no significant contribution of the res- 
onance decay 0~ — > E,*® 30 tt~ and set a combined upper 
limit B[Q-(Q+) -> 3*530(3*530)^] < 7.0 x 10~ 5 at 90% 
C.L. 
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